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Using the Density Matrix Renormalization Group and two-leg ladders, we investigate an electronic 
two-orbital Hubbard model including plaquette diagonal hopping amplitudes. Our goal is to search 
for regimes where charges added to the undoped state form pairs, presumably a precursor of a 
superconducting state. For the electronic density p = 2, i.e. the undoped limit, our investigations 
show a robust (7r,0) antiferromagnetic ground state, as in previous investigations. Doping away 
from p = 2 and for large values of the Hund coupling J, a ferromagnetic region is found to be stable. 
Moreover, when the interorbital on-site Hubbard repulsion is smaller than the Hund coupling, i.e. 
for U' < J in the standard notation of multiorbital Hubbard models, our results indicate the 
coexistence of pairing tendencies and ferromagnetism close to p = 2. These results are compatible 
with previous investigations using one dimensional systems. Although further research is needed to 
clarify if the range of couplings used here is of relevance for real materials, such as superconducting 
heavy fermions or pnictides, our theoretical results address a possible mechanism for pairing that 
may be active in the presence of short-range ferromagnetic fluctuations. 



PACS numbers: 71.10.Fd,71.27.+a,74.20.-z 
I. INTRODUCTION 

It is widely believed that magnetism is a fundamen- 
tal ingredient to explain the origin of high-temperature 
superconductivity in several materials. In fact, there 
is experimental evidence that the superconductivity in 
many heavy fermion (HF) compounds is mediated by 
spin fluctuationsii"— Mechanisms for superconductivity 
based on antiferromagnetism have been extensively dis- 
cussed for the Cu-based high-temperature superconduc- 
tors as wellj^ Recently, considerable excitement has been 
generated by the discovery of high-temperature supercon- 
ductivity in the iron pnictides.— Except for the cuprates, 
the iron-based superconductors now have the highest su- 
perconducting (SC) critical temperature Tc of any ma- 
terial (see for example Ref. Q). As in HF systems and 
cuprate superconductors, in the pnictides there is also 
evidence that the superconductivity is not mediated by 
the electron-phonon interactioui^i^ 

Magnetism and superconductivity can appear in dif- 
ferent ways. In some HF compounds, superconductiv- 
ity and antiferromagnetic (AFM) order co-exist?i while 
for the cuprates the superconductivity emerges after the 
long-range AFM order is destroyed by doping)^ In some 
HF systems, it is the superconductivity and ferromag- 
netism (FM) (as opposed to AFM order) that co-exist 
In this work, we will be interested in detecting a clear ev- 
idence of pairing of extra charges that are added to the 
undoped limit where magnetic order exists. Using a two- 
orbital model and two-leg ladders, it will be shown that a 



possible region for the robust coexistence of (spin triplet) 
pairing together with magnetic order is where ferromag- 
netism develops. This is in qualitative agreement with 
previous investigations carried out using one-dimensional 
systemsiii Our effort should be considered simply as pro- 
viding the first steps in relating pairing and magnetism 
in a complex two-orbital model via computational tech- 
niques on ladder geometries. Antiferromagnetic order, as 
found in the pnictides, could also be favorable for pair- 
ing tendencies in the spin singlet channel, as discussed 
recentlyJ^ii^ 

In principle, a theoretical investigation based on model 
Hamiltonians for strongly correlated materials starts 
with an effective tight-binding model, containing the 
minimum ingredients to describe the physics of the mate- 
rials under investigation. However, even if a well-defined 
reasonable model is used, it is still highly non-trivial to 
extract the ground state properties of these effective mod- 
els in two or three dimensions using unbiased numerical 
methods. In fact, at present there are no accurate tech- 
niques to study Hubbard-like models in dimensions two 
and three. Thus, in order to get at least some insight 
on the ground states properties of theses models, it is 
common practice to study the model Hamiltonians in 
quasi-one-dimensional geometries. In particular, a very 
popular route that has been used for several theoretical 
investigations is to study strongly correlated electronic 
systems using "ladder" geometries.— The A^-leg ladders 
consist of N chains of length L coupled by some param- 
eter (as, for example, fermionic hopping terms). The 
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two-dimensional system can in principle be obtained by 
considering the limits of both N and L sent to infinity, 
although in practice this is difficult to do. This ladder- 
based procedure has been used to investigate models for 
the high temperature superconductors^^ and for the HF 
systems4^ Some important results were obtained with 
this method. For example, research based on microscopic 
models for the high Tc superconductors^ as well as re- 
search on HF models^ indicate that superconductivity 
mediated by antiferromagnetic fluctuations can be stabi- 
lized, in agreement with several experiments. Thus, the 
use of ladders appears to be an important ingredient to 
unveil dominant ground state tendencies. Moreover, the 
hopping amplitudes that will be used in our investiga- 
tions below include next-nearest-neighbor diagonal hop- 
pings that are only possible when plaquettes exist in the 
lattice under consideration. 

Note that microscopic models that may present su- 
perconductivity induced by antiferromagnetism, such as 
the one-orbital Hubbard model and the Kondo Lattice 
model, have been extensively studied by several au- 
thors. However, microscopic models for superconductiv- 
ity in a ferromagnetic spin background have been much 
less explored, with the exception of studies using one- 
dimensional chains This may be caused in part by the 
perception that superconductivity and ferromagnetism, 
as opposed to antiferromagnetism, can not coexisti^i 
However, this perception has been challenged by the dis- 
covery of superconductivity and FM in the HF com- 
pounds UGea^a and URhGeJ^ Moreover, SC and FM 
were also observedi^ in the d-band metal ZrZn2. 

Motivated by the discovery of superconductivity in the 
HF compound UGe2 , a few years ago Karchev et al. pro- 
posed a one-band model to study the coexistence of su- 
perconductivity and ferromagnetismi^S However, other 
researchers have argued that the treatment used to inves- 
tigate that model, as well as the model itself, were not ap- 
propriate to describe the coexistence of these phasesi^"— 

Due to the lack of studies of microscopic models for su- 
perconductivity mediated by ferromagnetic fluctuations 
beyond one-dimensional chainSfii in this work we have 
decided to investigate a microscopic ladder model where 
the coexistence of superconductivity and ferromagnetism 
appears possible. In fact, it will be shown below that 
the two-orbital model that has been originally proposed 
to describe the low-energy physics of the iron-based 
supcrconductorsi^ii^i^i actually leads to the coexistence 
of pairing and spin ferromagnetic tendencies. Qualita- 
tive our results are compatible with those reported using 
chainsiii Although the model considered here may not be 
a proper effective model for superconducting ferromag- 
nets such as UGe2, we believe that the mechanism that 
bounds together the charges carriers (see section HI.B) 
is so simple and generic that our calculations may also 
apply in a variety of other models as well. 



II. MODEL 

In these studies, we have considered the following 
Hamiltonian defined on a two-leg ladder geometry 

H = - E (4.7.,Arf,+i..v,v + H. c.) 

- E ^7,7' (4.7<^.1^J\7V,2 + H. C.) 

Pj,x,\Pj,y,\ 

j,7.-^ 

~ 2"^E ' ^J,y,^ 

J, A 

+ ^EKWA^.^'^i.A^i^^'i^Arfi^yt^A + ff-c), (1) 
i,A 

where ^ creates an electron with spin projection a 
in the orbital 7 = x,y (d^z and dyz, respectively) at 
the rung j and leg A = 1,2, Sj,ry^\ is the electron spin 

density operator, Pj.^a.x = 4,7<t,a^j,7o-,A: and pj.^^x = 

The hopping amplitudes arc: t^]^ = iy,y ~ — <i, ty'y = 
ix,x = —t2i tl^'^ = i^'Jy = — ^3, t^'y — — i^'^y = —^4, if 
7 7^ 7', and zero otherwise. To avoid a proliferation of ex- 
tra parameters in our analysis, we have decided to fix the 
values of these hoppings from considerations previously 
used in the pnictide context. Our use of models originally 
devised for pnictides is simply based on the pragmatic 
observation that some pairing tendencies and ferromag- 
netic regions at large J were already observed in those 
modelsji^iH It should be clear though, that our research 
is mainly motivated by heavy fermion phenomenology. 

Following this strategy, then the hopping amplitudes 
ti, t2, ^3, and t4 are obtained using the Slater-Koster 
tight-binding scheme, and they are given fa yi^iiSi^s 

ti = -2{b^ -a^ + g^) /Apd 

t2 = -2{b^ -a^ -g^) /Apd 

ts = -{b^ + a^~g^)/Apd 

ti = -{ab-g^) /Apd, (2) 

where the Fe-As hopping amplitudes are a ~ 
0.324(pdcr) - 0.374(pd7r), b = 0.324(pd(T) + 0.123(pd7r), 
and g ~ 0.263(pd(T) -I- 0.31 (pdTr). A^d is the energy 
difference between the p and d levels. We have set 
[pda)"^ / Apd = 1 to fix the energy scale, and we use 
pdn/pda = —0.2, as previously discussedJ^ii^ Regard- 
ing the couplings U, J, and J', note that they are not 
independent, but they are assumed to satisfy the rela- 
tion U = U' + 2 J, which is strictly valid within a cubic 
environment for the full t2g sector i^i^ For an explicit 
derivation of this relation in the case of manganites see 
Ref. M 
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We have investigated the model defined above using a 
two-leg ladder of size 2xL, by means of the Density Ma- 
trix Renormalization Group (DMRG) technique,— un- 
der open boundary conditions (OBC), and keeping up 
to m=1400 states per block in the final DMRG sweep. 
We have carried out ^6—13 sweeps, and the discarded 
weight was typically 10^^ — 10~^ at the final sweep. In 
our DMRG procedure, the center blocks are composed 
of 16 states. We have used a FORTRAN DMRG code 
to calculate most results, and a C-I-+ DMRG code for 
additional validation)^ We have also confirmed some of 
our results against Lanczos Exact Diagonalization tech- 
niques, when possible. 

In this work, we will focus on the region of Hubbard 
and Hund parameters where Ucs = U — 3J = U' — J < 0, 
although a few results for positive values of Ucs will 
be presented as well. As shown below, in the region 
where UcS < there is a robust evidence of binding of 
holes/electrons close to density p = 2. In this regime 
U' < J. This inequality may not seem realistic at first 
sight, since the on-site spin triplet formation favored by 
J (Hund's rules) has its origin in the alleviation of the 
Coulombic energy penalization caused by U' . However, 
while in the full five-orbital model for Fe-based com- 
pounds U' < J is unphysical for the reason stated above, 
the two-orbital model is an effective model and it is un- 
clear how the main parameters arc affected by the pro- 
jection from five to two orbitals. Thus, the regime U' 
comparable to J may not be unrealistic. Clearly addi- 
tional investigations are needed to analyze if this regime 
of couplings is of relevance for real materials, such as 
heavy fermions. Ab-initio calculations are needed for this 
purpose (beyond the scope of the present analysis). 



III. RESULTS 

A. Magnetic Properties at p = 2 

Let us start our analysis by focusing on the density p = 
2. To investigate the magnetic order at this density, the 
Fourier transform of the real-space spin-spin correlation 
function was measured: 
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5j_j,_A- In Fig. [TJa), the spin struc- 
ture factor S{q) of the two-leg model is presented for 
several system sizes. As can be observed in this figure, 
there is a robust peak at wavevector q=(7r, 0), showing 
the tendency toward a stripe-like AFM order. This is the 
analog of the magnetic order found in pnictides but us- 
ing a two-leg ladder geometry. Similar results, obtained 
with Exact Diagonalization on small clusters, were re- 
ported beforej ^^i^^ Note that the peak increases with the 
system size, suggesting the development of a true long- 
range magnetic order in systems with a higher dimension. 
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Figure 1: (Color online), (a) Spin-structure factor S{q^) vs. 

for the two-leg ladder system with sizes L = 12, 16, and 
24, and for density p = 2, J = 1.5, and C/off ~ —0.5. (b) 
S{n, 0) as a function of J for ladders with linear sizes L — 8 
and UcH ~ —0.5. The inset shows the magnitude of this peak 
as a function of Ucti, for the coupling J — 1.5. (c) Spin-spin 
correlation Cspin(j) vs. j along the long ladder direction, for 
a system with size L = 24, and for the couplings C/cff ~ —0.5 
and J = 1.5. The dashed curve is a fit given by Eq. ([5]). 



The results for the spin structure factor show that 
along the y(x)-axes the spins are aligned following a 
ferromagnetic (AFM) order, at least at short distances. 
This stripe-like AFM structure is present in a wide range 
of parameters, as shown in Fig. [Ub) (inset), including 
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UeS > 0. Neutron scattering measurements for pnic- 
tides also show a similar spin order We have observed 
that the stripe-like AFM structure appears only when 
the plaquette-diagonal hopping amplitudes {t^ and ti) 
are of value similar as those of the nearest-neighbor hop- 
ping amplitudes. If we force t3= t4=0 then the peak 
in the spin structure factor S'(q) appears at wavevcctor 
q=(7r,7r). Note also that for the two- leg geometry, the 
(tt, 0) AFM state is not, naturally, degenerate with the 
(0,7r) AFM state. Due to this fact, a study in a ladder 
geometry could make a better connection with the two- 
dimensional results of the pnictide materials where (tt, 0) 
is favored over (0, tt) by a lattice distortion. 

For quasi-one-dimensional systems, a true long-range 
magnetic order is replaced by a power-law decay of the 
spin-spin correlations. Thus, in order to analyze the 
range of the magnetic order, we have also investigated 
the spin-spin correlation function along one of the legs 
(say, leg 1) defined as 
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where M is the number of site pairs satisfying I — 
— In practice, we have averaged over all pairs of sites 
separated by distance I, in order to minimize boundary 
effects (a few sites at the edges were also discarded while 
implementing this averaging procedure). 

In Fig. [IJc), the spin correlation function Cspin(j) is 
shown for the 2x24 cluster and using Ucs = —0.5. The 
dashed line is a fit of the numerical data with the function 



cos(7ra;) 



(5) 



Similar results were found for the 2x16 cluster. The ob- 
served power-law decay suggests that a two-dimensional 
system with the same model and parameters would de- 
velop long-range magnetic order at zero temperature. 



B. Doping with two holes or electrons 

Let us now consider the effect of doping with charges 
this ladder system. If tendencies toward the pairing of 
the extra charges are unveiled, they would be an indicator 
that this model could become superconducting in a two 
dimensional geometry, for the couplings here considered. 

Let us start with the calculation of the binding energy 
of two doped holes/electrons. This binding energy is de- 
fined as Ab = E{2) + E{0) - 2£'(1), where E{n) is the 
ground state energy with (4L + n) holes/electrons (4L is 
the number of electrons corresponding to the "undopcd" 
limit where there is an electron per orbital and per site) . 
On a finite system the binding energy Ab > is posi- 
tive if the electrons/holes do not form a bound state,— 
while in the thermodynamic limit Ab should vanish in 
the absence of pairing. On the other hand, if the extra 
holes/electrons form a bound state, then Ab < even on 



a finite cluster, and this would be indicative that effective 
attractive forces are present in the system. 

Before presenting our numerical results for the binding 
energies, let us first consider the following strong coupling 
regime defined by —Ucs = 3J — U = J — U' ^ 1. As 
it will be argued below, and as it was discussed in pre- 
vious investigations for one-dimensional chains?^ in this 
regime the binding of hole/electrons is clearly present. 

For completeness, let us address this limit in detail, 
although it is clear that large J compared with U' leads 
to an effective attractive interactioni^i Let us denote the 
states of the one-site problem via the symbol (^7^) i where 

the "arrows" Sx{sy) represent the electrons (with their 
spins projections) at the orbital x{y). For this one-site 
problem with two electrons and the limit considered here 
where J is large, the ground state energy is degener- 
ate. Its value is 62 = Ucs, and the three corresponding 

eigenstates are (^^^, and [(|)+ (^)]/%/2. Since 

—C/off ^ 1, we can approximate the Hamiltonian as 
H = ilkinetic + ^^int ~ ^^int, and the (highly degener- 
ate) ground state energy of the two-leg model becomes 
E{0) = 2L X 62 = 2L X U^s. 

Below, the dominant spin arrangement in the ground 
state at density p = 2 is shown to guide the discussion 




using only the up and down projections of the spin one 
states at every site for simplicity. Such a state is to be 
expected since for J ^ 1 the alignment of the two spins 
at the two orbitals in the same site will occur, and as 
U ^ 2J ^ 1 then having two electrons in the same or- 
bital is not allowed. Of course, any other configuration 
(such as a fully ferromagnetic state) is also equally likely 
as the one shown in the figure for —Ucs ^ 1. However, 
the hopping terms will lift the large degeneracy, and our 
numerical results for the spin-spin correlations presented 
before indicate that hopping terms favor the (tt, 0) AFM 
configuration, for a wide range of couplings in the un- 
doped limit. For this reason, here we have chosen to 
present the dominant spin arrangement of the ground 
state as a stripe-like AFM state, but it could be fully 
FM as well. Regardless of this detail, the arguments we 
will use below to obtain the energies E{n) will hold true 
both for a stripe-like AFM state, as well as for a FM 
state. 

Let us now add two extra electrons to the undoped 
system. In the limit being considered here, where the 
hopping amplitudes are negligible, the best way for the 
system to minimize its energy is to have the two electrons 
located on the same site, since in this way less on-site 
ferromagnetic links are broken. The on-site Hubbard U 
energy penalization is the same whether the doubly oc- 
cupied orbitals are at the same site or not, and since U' is 
negligible with respect to J in the limit considered, then 
the (effectively attractive) Hund coupling determines the 
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location of the extra charge, leading to the double occu- 
pation of both orbitals at the same sitei^i In this case, 
the dominant spin arrangement of the ground state is 
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Once again, a stripe-like AMF background is used, since 
for J ^ 1 and close to the density p ~ 2 our numerical 
data show (see Fig. |4|a)) that the (7r,0) AFM order is 
the dominant one, but it could have been FM as well. 

In the limit of couplings considered here, the ground 
state energy for the doped two-electron system is 

i;(2) = (2L + 5)C/cff + 8J. (6) 

Note that if the two extra electrons that are in the state 
(tt) ^^^^ move in opposite directions (after consid- 
ering the presence of small but nonzero hopping terms) 
they would break two ferromagnetic on-site links. Such a 
state would have a large energy and it is therefore "forbid- 
den" for — C/off ^ 1- However, if the two extra electrons 
move "together" (i.e. j) in the same direction, forming a 
spin triplet, then no other on-site FM links are broken. 
Thus, to minimize the energy the two-electrons added to 
the system must form a bound state, at least in the limit 
where the hoppings amplitudes are very small compared 
with the Hubbard and Hund couplings (the same spin 
dominant picture works when instead of adding electrons 
we remove two electrons, i.e. for the two- hole problem). 
This argument, known from previous investigations using 
chains explains the binding of electrons in the limit 
— J/eff ^ 1- More explicitly, the binding energy of the 
doped two-electrons/holes system in the limit considered 
here is given by 



|C/cff|. 



(7) 



Based on the discussion above we conclude that for 
— C/off S> 1, there is an indication of pairing, and perhaps 
superconductivity, in the two-orbital model. In fact, the 
on-site interorbital pairing state found here in this ex- 
treme regime results to be a spin triplet and transforms 
according to the irreducible representation of the 
group Z?4h'^ Interestingly, Exact Diagonalization calcu- 
lations in two-dimensional clusters, still in the FM state 
of the model but U' > J, found indications of a pairing 
state with the same characteristics and symmetry but 
with electrons at distance of one lattice spacing from each 
otherii^ii^ The argument presented above for this pair- 
ing is actually valid in any dimension, and also it works 
for the two doped holes case, as already mentioned. To 
confirm these argumentations, Ab was calculated numer- 
ically for large values of —C/off. An excellent agreement 
between the numerical data and the analytic expression 
(Eq. (O) was found. 

Now the crucial question is whether there is binding of 
holes/electrons for values of —C/off that may be of more 
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Figure 2: (Color online). (a)-(b) Phase diagram (J vs. 
— UcB ~ 3J — U) for the 2x3 cluster showing the region where 
the binding energies of electrons (a) and holes (b) is smaller 
than Ab < —0.1, indicative of pairing. Above the red line, 
the total spin saturates to its maximum value. Below this 
line, the total spin changes continuously from the maximum 
value to zero at J—Q. The blue lines show the region of J 
where Ab < —0.1 for UcS = —0.5. These lines are the same 
as presented in Fig.|3l (c) Ab vs. 1/L for some couplings (see 
legend). 



relevance for real materials. The answer to this ques- 
tion appears to be positive. In fact, the binding of elec- 
trons/holes has been observed numerically for C/off < 
and several values of the coupling J, as shown in Fig. 2M- 
More specifically. Figs. 2(a)-(b) show the region in the 
J — (3J— C/) plane where the binding energies of electrons 
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Figure 3: (Color online), (a) Phase diagram of the two-leg lad- 
der model Hamiltonian defined in Eq. 1. The region above the 
symbols is a fully-saturated ferromagnetic phase (see text). 
The blue lines correspond to regions where Ab < —0.1 for the 
cluster 2x3 with two doped electrons/holes, as presented in 
Figs. 2(a)-(b). 



(Fig.[5][a)) and holes (Fig.[5][b)) are less than -0.1, for the 
case of a 2x3 cluster. The region where Ab < —0.1 was 
chosen to be represented, as opposed to Ab=0, since pre- 
vious experience in the context of the cuprates^ suggests 
that this procedure effectively takes into account size ef- 
fects better. In practice, other values for this "cutoff" do 
not alter our qualitative conclusions. Similar results were 
found also for the 2x2 cluster. 

Thus far, only small clusters have been considered be- 
cause the numerical analysis of large clusters would be 
too time-consuming, particularly with regards to calcu- 
lating the hundreds of points that are required to extract 
comprehensive phase diagrams. However, larger system 
sizes were considered for a few selected sets of couplings, 
as shown in Fig. [2jc). In this figure, Ab vs. 1/L for 
some couplings is presented. Here, it is clearly observed 
that in the bulk limit the binding energies of added elec- 
trons/holes converge to nonzero values for some coupling 
sets. Close to p = 2, these results strongly indicate that 
there are pairing tendencies for L^off ^ . Thus, to the 
extend that future investigations show that J comparable 
to U' is a realistic regime for effective two-orbital mod- 
els, this provides a possible mechanism for pairing in real 
materials. 

In Figs. [5l^a)-(b), the magnetic phase diagram for the 
case of two doped electrons/holes is also presented. The 
region above the red (bold) line is a ferromagnetic phase 
with the maximum total spin 5totai = 2L— 1. Below this 
line, we have observed that the total spin changes con- 
tinuously from zero, at J = 0, up to its maximum value 
at Jc- As observed in these figures and for a large region 
of couplings, pairing (and presumably superconductivity) 
co-exists with ferromagnetic tendencies. 



C. Phase Diagram 

In Fig.|3]the phase diagram (J vs. density) of the two- 
leg ladder for Uctf = — 0.5 and several system sizes (see 
legend) is presented. For p = 2, it was observed that the 
total spin is zero for L even, and that it can be 0, 1, or 2 
for L odd, depending of the values of J. Note that for L 
odd and p = 2, the (7r,0) AFM configuration (H — I — h-) 
does not have the same number of + and - spins. 

For 2 < /9 < 2.5 (1 < p < 2), a ferromagnetic phase 
was found (the region above the symbols) with magnetic 
moment per site given by to = 2 — p/2 (to = p/2). The 
symbols indicate the value of Jc where the total spin sat- 
urates. The critical value Jc was determined by the level 
crossing of the energies in the sector with — •S'^j^j and 
= 5'[2tai ~ 1- Using this procedure, we were able to 
obtain Jc for large systems. For densities in the ranges 
p < 1.25 and p > 2.5, we have not found any trace of 
ferromagnetism. Below the FM region, it is very hard 
numerically to determine the total spin with good accu- 
racy for large systems. However, our results for the 2x2 
and 2x3 clusters with two doped electrons/ holes suggest 
that the total spin changes continuously, from maximum 
value at Jc to zero at J=0. The total spin can be ex- 
tracted from the spin structure factor at q = (0,0). For a 
few sets of couplings we also observed, through the value 
of S{q = (0, 0)), that in fact the total spin of the ground 
state changes continuously for the 2x8 cluster as well. 
These results suggest that the total spin varies continu- 
ously below the FM region present in Fig. [3] for any clus- 
ter sizes. Overall, our results are qualitatively compati- 
ble with those found in one-dimensional systemsJ^ Note 
also that tendencies to FM states at robust J were also 
reported via Exact Diagonalization methods on small 
clusters ji^ii^ 

We believe the FM phase is stabilized by a mechanism 
that has the same characteristics as the Double Exchange 
(DE) mechanism)^ In the original DE scenario, there 
are mobile and localized degrees of freedom. In the DE 
mechanism these degrees of freedom are separated and 
well defined. Although we do not have localized degrees 
of freedom in our model, from the perspective of one 
electron at a given orbital an electron at the same site but 
the other orbital behaves in some respects as a localized 
spin. For doped systems, when an electron moves from 
one site to the other, in order to minimize the kinetic 
energy and the energy related with the Hund coupling, 
all spins have to be aligned. 

The blue lines in Fig. [3] are the same that were pre- 
sented in Figs. [2Ia)-(b). We expect that the region of 
binding extends beyond these lines up to the density 
p — 2, forming regions in parameter space where super- 
conductivity exists inside the phase diagram. 

We have also measured the spin structure factor S'(q) 
away from the undoped density p = 2. In Fig. |4l S'(q) 
is presented for some particular densities for the two-leg 
ladder model with size L = 16, J = 1.5, and UcS = —0.5. 
As can be observed in Fig. IDJa), there is still a peak at 
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Figure 4: (Color online). Spin-structure factor S{qx) vs. g^; 
for the two-leg ladder system with size L = 16, J = 1.5, and 
f/cff = —0.5. (a) S{qx) for the densities p = 15/8, p = 33/16, 
and p = 17/8 (see legend), (b) S{qx) for the densities p = 7/4 
and p = 9/4. 

q = (ti", 0) for densities close to p = 2. Note that these 
peaks have smaller intensity than those found for p ~ 2 
in Fig. HJa), for the system with size L — 16. We have 
also observed that the height of the peak at q = (tt, 0) 
increases with the system sizes for the densities close to 
p = 2. These results indicate that a stripe-like AFM 
magnetic order also exists for densities close to p = 2. 
As shown in Fig. |31Jb), this order does not exist anymore 
for p > 2.2 and p < 1.7, at least within the precision 
of our calculations, and it is replaced by ferromagnetic 



tendencies. Note that for the electron doped case, there 
is a small peak at q = (0, tt) for densities p > 2.2. 



IV. CONCLUSION 

Using ladders, we have studied analytically and nu- 
merically a two-orbital Hubbard model. Via the DMRG 
technique we were able to investigate the model defined 
on a two-leg ladder geometry for systems with linear sizes 
up to L = 24. Our spin structure factor data show that 
for the "undoped" density p = 2, a stripe- like AFM order 
is present, as observed in previous Exact Diagonalization 
studiesji^ii^ We have also presented evidence for triplet 
pairing tendencies of added electrons/holes close to the 
density p = 2, in some range of couplings, in qualitative 
agreement with previous investigations using chains^ 
and with Exact Diagonalization calculations in a less ex- 
treme FM regime of models for pnictidesJ^ii^ More pre- 
cisely, we have found that pairing (and presumably su- 
perconductivity) and ferromagnetism co-exist for a large 
region of parameters in the regime U' < J. Even for U' 
comparable to J our results still indicate a (mild) ten- 
dency to pairing. Whether this range of couplings for U' 
and J is realized in real materials, such as heavy fermions 
or pnictides, is a matter to be decided via experiments, 
or with the help of ab-initio computer simulations. 
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